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Abstract 

We consider a gyroscopic system with two degrees of freedom under the action of small dissipative and 
non-conservative positional forces, which has its origin in the models of rotating bodies of revolution being 
in frictional contact. The spectrum of the unperturbed gyroscopic system forms a "spectral mesh" in the 
plane "frequency - gyroscopic parameter" with double semi-simple purely imaginary eigenvalues at zero 
value of the gyroscopic parameter. It is shown that dissipative forces lead to the splitting of the semi-simple 
eigenvalue with the creation of the so-called "bubble of instability" - a ring in the three-dimensional space 
of the gyroscopic parameter and real and imaginary parts of eigenvalues, which corresponds to complex 
eigenvalues. In case of full dissipation with a positive-definite damping matrix the eigenvalues of the ring 
have negative real parts making the bubble a latent source of instability because it can "emerge" to the 
region of eigenvalues with positive real parts due to action of both indefinite damping and non-conservative 
positional forces. In the paper, the instability mechanism is analytically described with the use of the 
perturbation theory of multiple eigenvalues. Explicit conditions are established for the origination of the 
bubble of instability and its transition from the latent to active phase, clarifying the key role of indefinite 
damping and non-conservative positional forces in the development and localization of the subcritical flutter 
instability. As an example stability of a rotating circular string constrained by a stationary load system is 
studied in detail. The theory developed seems to give a first clear explanation of the mechanism of self- 
excited vibrations in the rotating structures in frictional contact, that is responsible for such well-known 
phenomena of acoustics of friction as the squealing disc brake and the singing wine glass. 

Keywords: gyroscopic system, dissipative and non- conservative perturbations, indefinite damping, veering and 
merging of modes, spectral mesh, friction-induced oscillations, disc brake squeal, acoustics of friction, singing 
wine glass 

1 Introduction 

An axially symmetric shell, like a wine glass can easily produce sound when a wet finger is rubbed around 
its rim or wall, as it was observed already in 1638 by Galileo Galilei in his Dialogues Concerning Two New 
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Sciences [JJ [5J 131] • This principle is used in playing the glass harmonica invented by Benjamin Franklin in 
1757, which he called "armonica", where, in order to produce sound, one should touch by a moist finger an 
edge of a glass bowl rotating around its axis of symmetry [TS1 [531 GUI [S3] ■ This remarkable phenomenon has 
been studied experimentally; see e.g. [40] . However, an adequate analytical theory for its description seems 
to be still missing. Another closely related example of acoustics of friction is the squealing disc brake [22l [26] . 
This mechanical system produces sound due to transverse vibrations of a rotating annular plate caused by its 
interaction with the brake pads. Despite intensive experimental and theoretical study, the problem of predicting 
and controlling the squeal remains an important issue gl [TUl [^21 [Ml [111 [311 [171 HSJ EJJ . Significant but still 
poorly understood phenomena are squealing and barring of calender rolls in paper mills causing an intensive 
noise and reducing the quality of the paper [35] . 

The presence of multiple eigenvalues in the spectra of free vibrations of axially symmetric shells and plates 
is well-known. Already Rayleigh, studying the acoustics of bells, recognized that, if the symmetry of a bell were 
complete, the nodal meridians of a transverse vibration mode would have no fixed position but would travel 
freely around the bell, as do those in a wine glass driven by the moistened finger |22j . This is a reflection of 
the fact that spectra of free vibrations of a bell, a wine glass, an annular plate and other bodies of revolution 
contain double purely imaginary semi-simple eigenvalues with two linearly independent eigenvectors. 

Rotation causes the double eigenvalues of an axially symmetric structure to split [2], The newborn pair 
of simple eigenvalues corresponds to the forward and backward travelling waves, which propagate along the 
circumferential direction [21 [H [9l [TBI [T9l EQ] . Viewed from the rotational frame, the frequency of the forward 
travelling wave appears to decrease and that of the backward travelling wave appears to increase, as the spin 
increases. Due to this fact, double eigenvalues originate again at non-zero angular velocities, forming the 
nodes of the spectral mesh in the plane 'eigenfrequency' versus 'angular velocity'. The spectral meshes are 
characteristic for example for the rotating circular strings, rings, discs, and cylindrical and hemispherical shells. 
The phenomenon is apparent also in hydrodynamics, in the problem of stability of a vortex tube [27] and in 
magnetohydrodynamics (MHD) in the problem of instability of the spherically symmetric MHD a 2 -dynamo 
03J. 

It is known that striking the wine glass excites a number of modes, but rubbing the rim with a finger or 
bowing it radially with a violin bow generally excites a single mode |15) . The same is true for the squealing 
disc brake [26l [44J [4SJ [46l [48l [ST]. For this reason, we formulate the main problem of acoustics of friction of 
rotating elastic bodies of revolution as the description of the mechanism of activating a particular mode of the 
continuum by its contact with an external body. 

In case of the disc brake, the frictional contact of the brake pads with the rotating disc introduces dissipative 
and non-conservative positional forces into the system [481 151j . Since the nodes of the spectral mesh correspond 
to the double eigenvalues, they are most sensitive to perturbations, especially to those breaking the symmetries 
of the system. Consequently, the instability will most likely occur at the angular velocities close to that of the 
nodes of the spectral mesh and the unstable modes of the perturbed system will have the frequencies close to 
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Figure 1: The spectral mesh of system ([TJ when S = v = n = 0. 
that of the double eigenvalues at the nodes. This picture qualitatively agrees with the existing experimental 

data mm mm snug. 

In the following with the use of the perturbation theory of multiple eigenvalues [33J (Ml [Ml E3] we will 
show that independently on the definiteness of the damping matrix, there exist combinations of dissipative 
and non-conservative positional forces causing the flutter instability in the vicinity of the nodes of the spectral 
mesh for the angular velocities from the subcritical range. We will show that zero and negative eigenvalues in 
the spectrum of the damping matrix encourage the development of the localized subcritical flutter instability 
while zero eigenvalues in the matrix of non-conservative positional forces suppress it. Explicit expressions 
describing the movements of eigenvalues due to change of the system parameters will be obtained. Conditions 
will be derived for the eigenvalues to move to the right part of the complex plane. Approximations of the 
domain of asymptotic stability will be found and singularities of the stability boundary responsible for the 
development of instability will be described and classified. The methodology developed for the study of the 
two-dimensional system will finally be employed to the detailed investigation of the stability of a rotating circular 
string constrained by a stationary load system. 

2 The spectral mesh of a two-dimensional gyroscopic system 

Consider an autonomous non-conservative system described by a linear differential equation of second order 

x + (2ftG + <5D)x + ((J3 2 - n 2 )l + kK + i/N)x = 0, (1) 

where a dot over a symbol denotes time differentiation, x S M 2 , and I is the identity matrix. The real 
matrices D = D T , G = — G T , K = K T , and N = — N T are related to dissipative (damping), gyroscopic, 
potential, and non-conservative positional (circulatory) forces with magnitudes controlled by scaling factors <5, 
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fl, k, and v respectively; (3 > is the frequency of free vibrations of the potential system corresponding to 
5 = Q = K = h , = 0. The parameters and variables of the system are assumed to be non-dimensional quantities. 
Without loss in generality we assume det G = det N = 1 . Equation |T]) appears as a two-mode approximation 
of the models of rotating bodies of revolution in frictional contact after their linearization and discretization 

Separating time by setting x(t) = uexp(Ai) we arrive at the eigenvalue problem 

Lu = 0, L = IA 2 + (2f7G + <SD)A + (/3 2 - + kK + i<N. (2) 

Applying the Leverrier-Barnett algorithm [TT] to the matrix polynomial L we find the characteristic polynomial 
of the system 

P(A) = A 4 + <5trDA 3 + (2(/3 2 + f7 2 ) + 5 2 dctD + KtrK)A 2 
+ (4Qp + S(/3 2 - fi 2 )trD + <5ft(trKtrD - trKD))A 

+ K 2 detK + KtrK(/3 2 -ft 2 ) + (/3 2 -ft 2 ) 2 + i/ 2 . (3) 

In the absence of dissipative, external potential, and non-conservative positional forces (5 = k = v = 0) the 
characteristic polynomial ([3]) corresponding to the operator Lo(il) = IA 2 + 2A57G + (0 1 — fi 2 )I, which belongs 
to the class of matrix polynomials considered, e.g., in [21) . has four purely imaginary roots 

a+ = i(3 + in, a+ = -%p + in 

\- = \- = -i/3-m. (4) 

In the plane (f2, ImA) these roots considered as functions of O originate a collection of straight lines intersecting 
with each other (Fig. []} - the spectral mesh 43J . Two nodes of the mesh at f2 = correspond to the double 
semi-simple eigenvalues A = ±i(3. At the other two nodes at f2 = ±f2<2 there exist double semi-simple eigenvalues 
A = 0. The range |f2| < f2<j = j3 is called subcritical for the gyroscopic parameter f2 since the zero eigenvalue 
corresponds to the divergence boundary 21J. 

The double semi-simple eigenvalue i/3 at fl = O = has two linearly-independent eigenvectors U! and u 2 

u 2 = 4^ ( 1 I • (5) 
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Any linear combination of these vectors is an eigenvector too. The eigenvectors are orthogonal u^Ui = u 2 = 
and satisfy the normalization condition uj Ui = u 2 n u 2 = (2/3) -1 . 

Under perturbation of the gyroscopic parameter fl = fig + Af2, the double eigenvalue i(3 into two simple 
ones bifurcates. The asymptotic formula for the perturbed eigenvalues is [38] 



\± -a i -AO^ 11 + f 22 _1_ ■ AO /(Al _ /22) 2 , f ( (c\ 

\p=ifj + iAil ±tAild + /12/21 (6) 



where the quantities /j,- are 



n=0,A=i/3 
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Since G is a skew-symmetric matrix, we find 



/n = / 22 = 0, /i2 = -/ai = », (8) 
and Ay = i/3 ± iAfi = i/3 ± ifi in agreement with ([4]). 

3 Conservative, dissipative, and non-conservative deformation of 
the spectral mesh 

Following the approach of [55], we consider a perturbation of the gyroscopic system L (£l)+AL(f2) by dissipative, 
potential, and non-conservative positional forces breaking the symmetries of the initial system. Assuming that 
the size of the perturbation AL(f2) = 5AD + kK + j/N ~ e is small, where e = ||AL(0)|| is the Frobenius norm of 
the perturbation at £1 = corresponding to the double eigenvalue i(3, the behavior of the perturbed eigenvalues 
for small £1 and small e is described by the following asymptotic formula [38] 

,± -a . n fll + /22 . .£ll+£22 



, • /(^(/ll - hi) + e ll - £22) 2 un , v , Q s 

± «Y ^ 1" ("J12 + ei2)(S2/ 2 l + e 2 l ), (9) 

where the coefficients /y are given by the expressions jJJ and ey are small complex numbers of order s 

eij = ujAL(0)u. i = 0SujDui + kuJKu; + i/ujNu. t . (10) 

With the vectors (O we find 

en = 4^ + ^22, £12 = ifrfia + ^12 + ^ 

£22 = 4 dll + p* n ' £21 = 4 dl2 + i fcl2 ~^' (n) 
where dy and fcy denote the entries of the matrices D and K. Substituting the quantities {5J and (fTTj) into 
equation ([9]) we obtain 

A* = a/? + 1 ^ k <y±Vc, (12) 

where 

ft-M2\ r 2 (p\—p2\ 2 , ^„o , ^ \ .j. 2trKD — trKtrD 



and /ii, /X2 and pi, P2 are eigenvalues of the matrices D and K, respectively, and thus satisfy the equations 

/i 2 - /xtrD + dct D = 0, 

p 2 - ptrK + dct K = 0. (14) 

Separation of real and imaginary parts in equation (|12J) yields the relations 

( trD \ 4 / trD\ 2 , , 2 

4(ReA+— 5 J -4(RcA+— Sj Rcc-(Imc) 2 = 0, 

4llmA-/5-— k\ +4 ( ImA - /3 - — Rec-(Imc) 2 = 0, (15) 



where 



From the equations (|15p explicit expressions follow for the real and imaginary parts of the eigenvalues 
originated after the splitting of the double semi-simple eigenvalue i/3 



ReA = + Mg 5± Rec + V(Rec) 2 + (Imc) 2 



ImA = /3+ ^ K± ^/- RCC+ V /(R 2 Cc)2 + (Imc)2 . (17) 

The formulas (fl2l) - (fl7| describe splitting of the double eigenvalues at the nodes of the spectral mesh due 
to variation of parameters, including those corresponding to dissipative and non-conservative positional forces. 
Such splitting obviously leads to the deformation of the mesh, in particular, to the veering and merging of 
eigenvalue branches. Although the veering phenomenon in the systems with gyroscopic coupling was studied 
both numerically and analytically, e.g in [H [5J EH EE1 [23 , the explicit expressions (fT^|) - (fT7)) for the splitting 
of the double eigenvalues due to action of forces of all types were not previously derived. Note that the approach 
used in our paper differs from that of the cited works because it is based on the perturbation theory of multiple 
eigenvalues [33l [38 , 39J . The advantage of our approach is the description of the spectrum of the perturbed 
system by means of only the derivatives of the operator with respect to parameters and the eigenvectors of the 
multiple eigenvalue calculated directly at a node of the spectral mesh. 

3.1 Conservative deformation of the spectral mesh 

We first study the influence of a conservative perturbation with the matrix K on the spectral mesh of the 
gyroscopic system. Substitution of 5 = v = into the formulas ([15]) and (|TG|) yields 



j m X-0-r ± ^± K j _^ = ^±_^±] ^ ReA = 0. (IN) 



When k 7^ 0, equation (|T8[) describes a hyperbola in the plane ImA versus fl. Location of the branches of the 
hyperbola with respect to the lines of the unperturbed spectral mesh ImA = /3 ± fl depends on the definiteness 
of the matrix K, see Fig. [2J Indeed, according to equation (fTS)) the hyperbola has the asymptotes 

hn\ = p+ Pl * p2 K±a (19) 

The asymptotes cross each other above the node (0, 0) of the non-deformed spectral mesh for trK > 0, exactly 
at the node for p\ = —p2, and below the node for trK < 0. The branches of the hyperbola intersect the axis 
51 = at the points 

f3l=P+ % K > P2 = P+ % K - (20) 
If the eigenvalues p\ t 2 have the same sign, the intersection points are always located above or below the node 
of the spectral mesh for K > or K < 0, respectively, see Fig. [Ha). In case when one of the eigenvalues pi,2 is 
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Figure 2: Conservative deformation of the spectral mesh (k > 0): K is positive-definite (a); K is positive 
semi-definite (b); K is indefinite (c). 

zero implying semi-dcfinitcness of the matrix K, one of the branches of the hyperbola goes through the node of 
the spectral mesh and the other crosses the axis S7 = above the node, if K is positive semi-definite, or below 
it, if K is negative semi-definite, Fig.[2jb). If K is indefinite, one of the intersection points /?i.2 is located above 
the node and another one below the node, as indicated in Fig. (He). 

Therefore, the conservative deformation of the spectral mesh does not shift the eigenvalues from the imagi- 
nary axis near the nodes (0, ±/3) and thus preserves the marginal stability. However, the deformation pattern 
depends on the definiteness of the perturbation matrix K. In particular, in the degenerate case when det K = 0, 
one of the eigenvalue branches, originated after the splitting of the double eigenvalue, always passes through 
the point corresponding to the node of the unperturbed spectral mesh, Fig. [3]^b). 

3.2 Creation and activation of the latent sources of instability by dissipation 

Now we consider the effect of dissipative forces on the stability of the gyroscopic system and study its dependence 
on the properties of the matrix D. Assuming v = k = in expression (j 13(1 we rewrite formula (|12j) in the form 

X = i0— ^±^5 ± 2 ~W . (21) 

Since Imc — 0, equations (|15p can be transformed into 

( ReA + f±±J» S ) 2 + O 2 = (Ml ~ M2)2 S 2 , lm\ = f3 (22) 
\ 4 ) 16 

when 

n 2 (/il-/i2) V < 0, (23) 
16 

and into 

Q 2 -(ImA-/?) 2 = (/il 7 M2) % 2 , ReA = -^±^S, (24) 
16 4 

when the sign in inequality (|23[) is opposite. For a given 8 equation (|24[) defines a hyperbola in the plane 
(f2,ImA), while ([22|) is the equation of a circle in the plane (fi,ReA), as shown in Fig. [Ha,c). For tracking the 
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Figure 3: Origination of a latent source of the subcritical nutter instability in presence of full dissipation: 
Submerged bubble of instability (a); coalescence of eigenvalues in the complex plane at two exceptional points 
(b); hyperbolic trajectories of imaginary parts (c). 

complex eigenvalues due to change of the gyroscopic parameter fi, it is convenient to consider the eigenvalue 
branches in the three-dimensional space (f2, ImA, ReA). In this space the circle belongs to the plane ImA = (3 
and the hyperbola lies in the plane ReA = —S(fJ-i + A t 2)/4, see Fig. |4ja,c). Below we show that the circle of the 
complex eigenvalues - the bubble of instability |23] - plays the crucial role in the development and localization of 
the subcritical flutter instability and that its properties depend on whether the matrix D is definite or indefinite. 

3.2.1 Full dissipation and pervasive damping: A latent state of the bubble of instability 



From the formulas (|22[) and (|24[) we see that the basic geometric characteristics of the trajectories of eigenvalues 
of the dissipatively perturbed gyroscopic system are defined by the eigenvalues /xi and p,i of the matrix D. For 
example, the radius of the bubble of instability r b and the distance d b of its center from the plane ReA = are 

Km - H2)5\ +^ 2 )(5| 
r b = , d b = . (25) 

As a consequence, the bubble does not intersect the plane ReA = under the condition d b > r&, which is 
equivalent to the inequality 

^ lAi2 = detD > 0. (26) 
Moreover, the bubble of instability is "submerged" under the surface ReA = in the space (SI, ImA, ReA), if 

<5(Mi+M2) = <5trD>0. (27) 

Since the inequalities (|26p and (|27[) imply positive semi-definiteness of the matrix 8D, we conclude that 
the role of full dissipation or pervasive damping is to deform the spectral mesh in such a way that the double 



Figure 4: The mechanism of subcritical flutter instability (bold lines): The ring (bubble) of complex eigenvalues 
submerged under the surface ReA = due to action of dissipation with det D > - a latent source of instability 

(a) ; repulsion of eigenvalue branches of the spectral mesh due to action of non-conservative positional forces 

(b) ; emersion of the bubble of instability due to indefinite damping with det D < (c); collapse of the bubble of 
instability and immersion and emersion of its parts due to combined action of dissipative and non-conservative 
positional forces (d). 
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semi-simple eigenvalue is inflated to the bubble of complex eigenvalues (|22p connected with the two branches of 
the hyperbola (jiM)) at the points 



and to plunge all the eigenvalue curves into the region ReA < 0. The eigenvalues at the points (|28[) are double 
and have a Jordan chain of generalized eigenvectors of order 2. In the complex plane the eigenvalues of the 
perturbed system move with the variation of Q along the lines ReA — —db until they meet at the points (|28[) 
and then split in the orthogonal direction; however, they never cross the imaginary axis, see Fig. [3jb). Note 
that a similar process of unfolding the semi-simple eigenvalue (diabolical point) into two double eigenvalues 
with the Jordan blocks (exceptional points) for Hermitian matrices under arbitrary complex perturbation and 
its application to crystal optics have been described in [501 13H] ■ 

The bubble of instability has two remarkable properties important for the explanation of the phenomenon of 
squeal. First, there can exist perturbations causing its growth and emersion above the surface ReA = (flutter); 
second, for small 6 the instability is localized in the vicinity of the frequency ImA = [3 and the value of the 
gyroscopic parameter O = 0. The narrow frequency band is a characteristic property of squeal. We conclude 
that the bubble is a latent source of subcritical flutter instability localized in a narrow range of change of the 
gyroscopic parameter with |f2| < U,d at the frequency corresponding to the double eigenvalue of the non-rotating 
system. 

3.2.2 Indefinite damping: An active state of the bubble of instability 

As it is seen from equations (|25|) , the radius of the bubble of instability is greater then the depth of its submersion 
under the surface ReA = only if the eigenvalues fj,± and [ii of the damping matrix have different signs, i.e. if the 
damping is indefinite. The damping with the indefinite matrix appears for example in the systems with frictional 
contact when the friction coefficient is decreasing with relative sliding velocity as well as in the problems of 
propagation of waves in dissipative media [17, 24, 26 . Indefinite damping is known to be a destabilizing factor 
[l8l [49l [50] . In our case it leads to the emersion of the bubble of instability meaning that the eigenvalues of the 
bubble have positive real parts in the range fl 2 < fl 2 r , where 



is found from the equation (|22|) after assuming there ReA = 0. The right side of the formula ([29]) is real only 
for detD < 0, i.e. for the indefinite matrix D. 

We see that the domain of asymptotic stability is defined by the constraints <5trD > and tt 2 > fl 2 r . In the 
plane of the parameters 6 and f2 for detD < it has a form shown in Fig. 02a). Due to the singularity existing 
at the origin, an unstable system with indefinite damping can be stabilized by sufficiently strong gyroscopic 
forces, as shown by the dashed line in Fig. [5][a). With the increase of detD the stability domain is getting 
wider and for detD > it is defined by the condition <5trD > 0, Fig. EKc). At detD = the line f2 — does 
not belong to the domain of asymptotic stability, Fig. 03[b). 



ImA = 0, ReA 



5(ni + vb)/4, Q = ±6(p 1 -fi 2 )/i, 



(28) 




(29) 
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Figure 5: Approximation based on the first-order perturbation theory of double eigenvalues to the domain of 
asymptotic stability (white) and its boundary (bold lines) for the dissipatively perturbed gyroscopic system (1) 
when v = k = 0, trD > and detD < (a); detD = (b); detD > (c). 

Therefore, changing the damping matrix <5D from positive definite to indefinite we trigger the state of the 
bubble of instability from latent (ReA < 0) to active (ReA > 0), see Fig. HJa,c). Since for small S we have 
£l cr < Qd, the flutter instability is localized in the vicinity of — 0. As we show below, the non-conservative 
positional forces play the similar role. 

3.3 Activation of the bubble of instability by non-conservative positional forces 

In the absence of dissipation non-conservative positional forces destroy the marginal stability of gyroscopic 
systems [6j [7]. One can easily check that this property is valid for system (1) by assuming 6 = n = in the 
formulas l|12p and ([T5)l . which yield 

A± = z/?±ztt±^, A± = -i/3±*n=Fp. (30) 

Equations l|30p show that the eigenvalues of the branches i/3 + i£l and —i/3 — iQ, of the spectral mesh 
get positive real parts due to perturbation by the non-conservative positional forces. The eigenvalues of the 
other two branches are shifted to the left from the imaginary axis, see Fig. 0Jb). Note that the real part of 
the perturbed eigenvalue is proportional to 1//3, which means that the destabilizing effect of non-conservative 
positional forces is less pronounced for higher frequencies. This is important in the study of systems with many 
degrees of freedom. Thus, the instability induced by the non-conservative forces only is not localized near the 
nodes of the spectral mesh, in contrast to the effect of indefinite damping. We show now that in combination 
with the dissipative forces, both definite and indefinite, the non-conservative forces can create subcritical flutter 
instability in the vicinity of diabolical points. 

Recall that the real and imaginary parts of the eigenvalues originated after the spitting of the double eigen- 
values due to combined action of dissipative and non-conservative positional forces are given by the expressions 
(fl5| and (fT7|) . where we assume k = Q. Multiplying equations (fT7l) we find that the trajectories of the eigenvalues 
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Figure 6: Subcritical nutter instability due to combined action of dissipative and non-conservative positional 
forces: Collapse and emersion of the bubble of instability (a); excursions of eigenvalues to the right side of the 
complex plane when Q goes from negative values to positive (b); crossing of imaginary parts (c). 



in the complex plane are described by the formula 



RcA + ^) (ImA-/3) = |. 



(31) 



When v = and 5 =/= is given, the eigenvalues move with the variation of the gyroscopic parameter 51 along 
the lines ReA = — trD/4 and ImA = (3 and merge at the points (|28|) . see Fig. [3{b). 

Non-conservative positional forces with v ^ destroy the merging of modes. As a consequence, the eigen- 
values move along the separated trajectories. According to relations (|30[) the eigenvalues with |ImA| increasing 
due to an increase in |f2| move closer to the imaginary axis then the others, as shown in Fig^b). Therefore, in 
the space (fi, ImA, ReA) the action of the non-conservative positional forces separates the bubble of instability 
and the adjacent hyperbolic eigenvalue branches into two non-intersecting curves, see FigUJd). It is remarkable 
that the form of each of the new eigenvalue curves carries the memory about the original bubble of instability. 
Due to this spectral peculiarity the real parts of the eigenvalues can be positive for the values of the gyroscopic 
parameter localized near = in the range fl 2 < fl 2 r . Taking into account that ReA = at the critical values 
of the gyroscopic parameter we find Q cr from the equations (|17|) 

trD / v 2 - S 2 f3 2 detD 



a 



-5 2 f3 2 (trB/2) 2 ' 



(32) 



Additionally, it follows from (1 1 7[) that the eigenfrequencies of the unstable modes from the interval Vt 2 < f2 2 r 
are localized near the frequency of the double semi-simple eigenvalue at the node of the undeformed spectral 
mesh: oj~ r < ui < oj+. 



aJ,. 



v 2 - S 2 (3 2 detD 



(33) 



2/3 Y *; 2 -£ 2 /3 2 (trD/2) 2 ' 

When the radicand in formulas (|32[) and (|33[) is real, the eigenvalues move in the complex plane making 
the excursion to its right side, as shown in Fig. [6^b). As it follows from (|32l) . in presence of non-conservative 
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positional forces such excursions behind the stability boundary are possible for the eigenvalues, even when 
dissipation is full (detD > 0). One can say that similarly to the indefinite damping the non-conservative 
positional forces activate the latent sources of flutter instability created by the full dissipation. 

The equation (|32[) describes the surface in the space of the parameters S, v, and f2, which is an approximation 
to the stability boundary separating the domains of asymptotic stability and flutter. For better understanding 
the shape of this surface we extract the parameter v in (|32[) 



If dot D > and f2 is fixed, the formula (CM)) describes two independent curves in the plane (<5, v) intersecting 
with each other at the origin along the straight lines given by the expression 



However, in case when damping is indefinite and det D < 0, the radical in the formula Q34p is real only for 
8 2 < — 4f7 2 / det D meaning that ([3^]) describes two branches of a closed loop in the plane of the parameters 5 
and v. The loop is smooth at every point except at the origin, where it is self-intersecting with the tangents 
given by the expression |35|) . Therefore, for a given O this curve looks like figure "8". When fi goes to zero, 
the size of the self-intersecting curve tends to zero too. We conclude that in case det D < the shape of the 
surface described by equation ([3"!?]) or (fM]l is a cone with the "8" -shaped loop in a cross-section, see Fig. [7{a). 
Due to self-intersections the cone consists of four pockets. The asymptotic stability domain is inside the two 
of them, selected by the inequality ^Ff\ . as shown in Fig. [7{a). The singularity of the stability domain at the 
origin is the degeneration of a more general configuration described first in |50j . 

The domain of asymptotic stability bifurcates when det D changes from negative to positive values. This 
process is shown in Fig. [Jj Indeed, according to formula (|34|) with increasing det D < two pockets of the 
domain of asymptotic stability move towards each other until they have a common tangent line v = at 
detD = 0, sec Fig. [TJb). When detD is positive, this temporarily glued configuration unfolds to the unique 
domain of asymptotic stability bounded by the two surfaces intersecting along the 17-axis, as shown in Fig. [7jc). 

In Fig. [7Ja) we see that in case of indefinite damping there always exists an instability gap due to the 
singularity at the origin. Starting in the flutter domain at SI — for any combination of the parameters 5 and v 
one can reach the domain of asymptotic stability at higher values of (gyroscopic stabilization), as shown in 
Fig. [7ja) by the dashed line. The gap is responsible for the subcritical flutter instability localized in the vicinity 
of the node of the spectral mesh of the unperturbed gyroscopic system. When detD = 0, the gap vanishes in 
the direction v — 0. In case of full dissipation (detD > 0) the singularity at the origin unfolds. However, the 
memory about it is preserved in the two instability gaps located in the folds of the stability boundary with the 
locally strong curvature, Fig. [7Jc). 

We see that in case of full dissipation for some combinations of the parameters S and v the system is 
asymptotically stable at any Q. There exist, however, the values of 5 and v for which one can penetrate the 




(34) 



v = ± 



/3trD 



(35) 
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a) 




Figure 7: Domains of asymptotic stability in the space (<5, is, f2) for different types of damping: Indefinite 
damping detD < (a); semi-definite (pervasive) damping detD = (b); full dissipation detD > (c). 

fold of the stability boundary with the change of il, as shown in Fig. [7Jc) by the dashed line. For such 8 and v 
the flutter instability is localized in the vicinity of f2 = 0. It is remarkable that in presence of non-conservative 
positional forces the system with full dissipation can suffer from the subcritical flutter localized near the nodes 
of the spectral mesh. A good illustration of this fact is the formula for the maximal real part of the unstable 
eigenvalue attained at f2 = (see Fig. [Ha)) 



From our previous considerations it follows that the phenomenon of the local subcritical flutter instability 
is controlled by the eigenvalues of the matrix D. When both of them are positive, the folds of the stability 
boundary are more pronounced if one of the eigenvalues is close to zero. If one of the eigenvalues is negative and 
the other is positive, the local subcritical flutter instability is possible for any combination of 5 and v including 
the case when the non-conservative positional forces are absent (y = 0). 

Even if the structure of the damping matrix D is unknown, we realize that the main role of dissipation of 
any kind is the creation of the bubble of instability. It is submerged below the surface ReA = in the space 
(Q, ImA, ReA) in case of full dissipation and partially lies in the domain ReA > when damping is indefinite. 
Non-conservative positional forces destroy the bubble of instability into two branches and shift one of them to 
the region of positive real parts even in case of full dissipation. Since the branch remembers the existence of 
the bubble, the instability is developing locally near the nodes of the spectral mesh. 

Therefore, the instability mechanism behind the squeal is the emersion (or activation) due to indefinite 
damping and non-conservative positional forces of the bubbles of instability created by the full dissipation in the 
vicinity of the nodes of the spectral mesh. 




(36) 
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Figure 8: A rotating circular string and its "keyboard" constituted by the nodes (marked by white and grey) 
of the spectral mesh (only 30 modes are shown). 

4 Example. A rotating circular string constrained by a stationary 
load system 

Following Yang and Hutton [16] , we consider a rotating circular string of displacement W(tp, r), constrained at 
if = by a stationary load system consisting of a spring, a damper, and a massless eyelet generating a constant 
frictional follower force F |19j on the string, as shown in Fig.[8l The parameters r and p are the radius and mass 
per unit length of the string. The rotational speed of the string is 7. The following assumptions are adopted in 
developing the governing equation of the problem: the circumferential tension P in the string is constant; the 
stiffness of the spring supporting the string is K and the damping coefficient of the viscous damper is D; the 
velocity of the string in the ip direction has constant value "fr [16] . 

Introducing the following non-dimensional variables and parameters 



t 




W 



P' 



Kr F 



d = 



D 



(37) 

(38) 
(39) 

The boundary conditions (I39|) reflect the continuity of the string displacement and the discontinuity of its slope 
following from the force balance at the eyelet in the assumption of smallncss of the norm of w and w v with 
respect to r. The inclusion in Fig. [8] shows the mutual configuration of the vectors of the frictional follower 
force F and of the circumferential tension P. 



we arrive at the non-dimensional governing equation and boundary conditions [16j : 

Wu + 2ftw tip - (1 - fl 2 )w vv = 0, 
w(0,t) - w(2-K,t) = 0, (1 - n 2 )[w v> (2Tr,t) - w v (0,t)] + kw(0,t) + dw t (0,t) + pw v (0,t) = 0. 
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Separation of time by the substitution w(ip, t) = u(ip) exp(Ai) yields the boundary value problem 

Lu = \ 2 u + 2n\u' - (1 - tt 2 )u" = 0, (40) 



MO! uVlirl = 0. »'(()! - ,/(2,t) = ^±£«(0) + Y^^'CO), (41 i 



where prime denotes partial differentiation with respect to ip. The non-self-adjoint boundary eigenvalue problem 
depends on four parameters O, d, k, and ji expressing the speed of rotation, and damping, stiffness, and friction 
coefficients of the eyelet, respectively. 

Taking the scalar product (Lu, v) — vLudtp, where the bar over a symbol denotes complex conjugate, 
then integrating it by parts and employing the boundary conditions (|41|) we arrive at the boundary value 
problem adjoint to ((3D]) and |gTJ 

L*v = \ 2 v - 2n\v' - (1 - fl 2 )v" = 0, (42) 

«(0) - v(2n) = - T ^ V (27T), «'(0) - v'(2n) = j^«(27r) + (1 ^^ ^tt). (43) 

Let us first consider the string without constraints (d = k = fx = 0). Then, the system is gyroscopic so that 
the eigenfunctions of the adjoint boundary value problems corresponding to a purely imaginary eigenvalue A 
coincide, i.e. v — u. Assuming the solution of the equation (|40[) in the form u = C\ exp (jz^i<pj +C2 exp ^ ^TH^) 
and substituting it into the boundary conditions (|41|) we obtain the characteristic equation 

-2ttAQ 

q 1 < ttA . ttA e " 2 -i 

8A sm sm ToT^^T = °' (44) 

The roots of the equation (|44p are 

A+ = m(l + fi), A"=m(l-n), (45) 

where n € Z. They are the eigenvalues of the boundary eigenvalue problem (|40|) , (|4ip with the eigenfunctions 

w+ = cos(n<p) — z sin(n^j), u~ = cos(n<p) + i sm(nip), (46) 

respectively. The eigenvalues are purely imaginary and form a mesh of lines intersecting with each other in the 
plane ImA versus J7, as shown in Fig. [5] 

The eigenvalues (|4"5j) are simple almost at all values of f2 excepting those corresponding to the nodes of the 
spectral mesh. Indeed, two eigenvalue branches A^ = in(l + eQ) and Af n = im(l + 5Q), where e, 8 = ±, intersect 



each other at £1 = Q eS 



inn 



= (47) 
mo — ne 



and originate the double eigenvalue 



s _ inm(5 - e) 



Knn — ——J 7" ) (48) 

mo — ne 



which has two linearly independent eigenfunctions 

u e n = cos(rvp) — ei sm(n<p), u s m — cos(m(p) — Si sm(mip). (49) 



1G 



The nodes (|4"T)) , (|4"8|l of the spectral mesh of the rotating circular string in the absence of the external 
loading are marked by white and grey dots in Fig. [5J At Q — the spectrum of the non-rotating circular string 
consists of the double semi-simple purely imaginary eigenvalues in, n £ Z, each of which splits into two simple 
purely imaginary eigenvalues due to change in the angular velocity [9J [10], [16l [20] . At £1 = ±1 all the eigenvalue 
branches cross the axis ImA = 0, see Fig. [8j In the following we will consider the spectrum for the angular 
velocities from the subcritical range when il £ (— 1, 1). 

Now we consider the deformation of the spectral mesh caused by the interaction of the rotating string with 
the external loading system. We proceed analogously to our investigation of a two-dimensional system and 
study the splitting of the double eigenvalues at the nodes of the spectral mesh. For this purpose we use the 
perturbation theory of multiple eigenvalues of non-self-adjoint differential operators developed in [551 1551 135) . 
According to this theory the perturbed eigenvalues are expressed by the following asymptotic formula 



A 



fee I f 



65 
mm 



f ee I f S5 
^nn ' ^mm 



± 



{fnn fH + e nn e 22) 



if, 



eS 



e eS )(fSe +£ 8e ) 
nm ' ^nmJ \J mn 1 ^mn / 



)(f s 

J\Jn 



c 6e 



(50) 



The coefficients f^ m are denned by 

xeS 

J nm 



An, 



while the quantities e^ m are 

*6 = ( dX nm + fcR(0)^(0) + ^< / (0)^(0) 

nm I — ~ ~ 

where Ail = il — il^m- Calculating the integrals in (f5T|) and ([52"]) and taking into account expressions 
we obtain 

fZ = -isnAil, f± = = 0, f™ m = -iSmAil, 



(51) 



(52) 
(g7J and 

(53) 



and 



dX 



eS 



k — sinfx 



dX 



eS 



k — einfi 



dX n s m + k — Simfj, 



d\„ m + k — Simfi 



mn a ■ i i mm a - 

VKl^Jnm VKTtll 

Taking into account the results of our calculations ([55]) and ([5JJ we find 



.sn + 5m 
A = X E JL + 1 An + 1 



n + m , £ s .n e + 6 



Sunm 



(dX% m + k) + 



8tt 



where 



/ i £n ~ Sm A r, , i m ~ n ^X eS | fc) | £ ~ 5 fJ \ ( dA »m + k ~ ien P){dX n 5 m + k - i8mp) 



8wmn 



8tt 



167r 2 nm 



(54) 



(55) 



(56) 
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According to the experimental data [35J 123 US ESI US] the frequency of sound emitted by a singing wine 
glass and a squealing laboratory brake at low spins is close to a double eigenfrequency of the non-rotating 
bodies. For this reason it is important to study first the influence of external stiffness, damping and friction on 
the deformation of the spectral mesh near the nodes corresponding to £1 = 0. 

At ft — the double eigenvalue X^m originates due to intersection of the eigenvalue branches A^ and X s ml 
where m = n and s = —5. Taking this into account, we find from (|55|) and (|56p the expressions describing 
splitting of the double eigenvalue in due to action of gyroscopic forces and an external spring 



k k 2 
X = in + i- — ±i\n 2 n 2 + -— — . (57) 
47m V lon z n z 

Equation ([57]) describes a hyperbola in the plane ImA versus f2, which intersects the axis f2 = at the values 
A = n and A = n + see Fig. Ufa). The gap between the branches decreases with the increase of the 
number n of a mode. It is remarkable that the lower branch passes through the point corresponding to the 
node of the spectral mesh of the non-perturbed gyroscopic system, which agrees with the numerical results of 
[16j . Remember that in case of two-dimensional systems the reason for such a degenerate behavior is a zero 
eigenvalue in the matrix K of external potential fores. 

Analogously, from expressions (|55p and (|56p follows the asymptotic formula for the eigenvalues originating 
after the splitting of the double eigenvalue in at f2 = due to perturbation by the gyroscopic forces and an 
external damper 

d j~d 2 

According to (|58|) the real parts of the eigenvalues as functions of f2 originate a bubble of instability in the plane 
ReA versus il 

j \ 2 ,2 

ReA H ) +n 2 n 2 = ImA = n. (59) 

in ) lon z 

The ellipse (|59l) is submerged under the plane ReA = in the space (f2, ImA, ReA) so that it touches the plane 
at the origin, as shown in Fig. Efb). The ellipse is connected with the branches of the hyperbola of complex 
eigenvalues 



\ = in- — ±\j Tr ^~n 2 VL 2 . (58) 



d 2 d 

n 2 n 2 - (ImA - nf = -, ReA = . (60) 

16n 2 4-7T 



The horizontal diameter of the ellipse decreases with the increase in the mode number n, while the vertical one 
does not change. As we see, the external damper creates a latent source of local subcritical flutter instability 
exactly as it happens in two dimensions when the matrix of dissipative forces D is semi-definite, i.e. when it has 
one zero eigenvalue. The range of change of the gyroscopic parameter corresponding to the latent instability is 
located compactly around the origin and decreases with the increase in n. 

The deformation of the spectral mesh near the double eigenvalue in at Q = due to combined action of 
gyroscopic forces and external friction is described by the expression 



\ = ,„±\ (inn+^Y-J^, (61) 
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Figure 9: Deformation of the spectral mesh of the rotating string near the nodes (0, 3), (0, 2), and (0, 1) caused 
by the action of the external spring with k = 0.3 (a), damper with d = 0.3 (b), and friction with /i = 0.3 (c), 
respectively. 
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following from (f53|) and (|56p. 

The imaginary parts of the eigenvalues of the deformed spectral mesh 

ImA = n ± -!- V 27r 2 n 2 fi 2 ± Trnfi V^n^ + £t 2 (62) 

27T V 

cross at the node (0, n), as in the non-perturbed case. However, the crossing is degenerate because the eigenvalue 
branches touch each other at the node, see Fig.[9jc). Indeed, expanding expression (|62|) in the vicinity of f2 = 
we find that 

ImA — n± — y^rn/Ltjfi] + 0{n 3/2 ). (63) 

27T 

Clearly, at = the imaginary parts do not split due to non-conservative perturbation from the eyelet. It 
is instructive to note that for Q — > oo the imaginary parts asymptotically tend to n(l ± fi). By this reason 
for small perturbations the spectral mesh depicted in the plane ImA versus SI looks non-deformed at the first 
glance. 

The real parts of the deformed spectral mesh described by the expression 

ReA = ± — J-2Tr 2 n 2 n 2 ± Trnn^/^n 2 ^ 2 + y? (64) 
2n 

cross at il = 0, so that ReA = 0. The crossing of the real parts is degenerate Fig. E^c), which is confirmed by 
the asymptotic expression 

ReA = ±^-^Trnn\n\ + 0{n 3/2 ). (65) 
For f2 — > oo the real parts follow the asymptotic law 

ReA = ± £ T m^ + °^ ^ 

As is seen in Fig. Etc), the real parts almost always are close to the lines ±/x/ (47r), except for the vicinity of the 
node of the spectral mesh, where the real parts rapidly tend to zero. This behavior agrees with the results of 
numerical calculations of [P6] . 

We see that the double semi-simple eigenvalue in does not split due to variation of only the parameter of non- 
conservative positional forces \x. In the two-dimensional case this would correspond to the degenerate matrix 
N, dctN = 0. Since N is skew-symmetric, the degeneracy means N = 0, i.e. the absence of non-conservative 
perturbation. In case of more then two degrees of freedom the degeneracy of the operator of non-conservative 
positional forces leads to the cuspidal deviation of the generic splitting picture registered in Fig. [9jc). 

We have considered local deformation of the spectral mesh of the rotating string near the nodes at f2 = 0. 
However, the formulae (|55[) and (|56p enable us to investigate similarly the deformation of the mesh in the 
vicinity of all the other nodes. 

Concluding, we note that the source of perturbation of the rotating string concentrated at one point leads 
to the deformations of its spectral mesh, which correspond to that caused by the semi-definite matrices of con- 
servative, dissipative and non-conservative forces in the two-dimensional case. Zero eigenvalues of the damping 
operator encourage the activation of the latent bubble of instability, while the zero eigenvalue of the operator of 
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non-conservative positional forces suppresses this process. To get local subcritical nutter instability, described 
in the previous sections for the finite-dimensional model, the operators of dissipative and non-conservative 
perturbations must be generic, which excludes their semi-definiteness. We notice that the singular sources of 
momentum and energy in the problem of stability of a rotating string as a reason for degeneracies were dis- 
cussed also in [28j [31] . One of the ways to avoid this degeneracy is to consider not pointwise but distributed 
contact with the dissipative, stiffness and friction characteristics depending on the material coordinates. A step 
in this direction is taken in |51j . where a model of distributed pads was developed. For simplicity, in |51j the 
characteristics of the pins constituting the pads were assumed not depending on the coordinates, which led to 
the same semi-definite degeneracy. This could be a reflection of the so-called Herrmann-Smith paradox of a 
beam resting on a Winkler- type elastic foundation and loaded by a follower force [35]. The degeneracy in the 
Herrmann-Smith problem is removed by assuming a non-uniform modulus of elasticity. Similar modification of 
the model of the distributed brake pads of [51] could give generic perturbation operators and open the way to 
the modeling of the disk brake squeal catching its most significant features. 

5 Conclusion: Disc brake as a musical instrument 

As we already mentioned, the principle of activating sound by friction is the same for a wine glass, a disc brake, 
and the glass harmonica. The latter is an ancient musical instrument for which the famous " Dance of the Sugar 
Plum Fairy" in the first edition of "The Nutcracker" ballet was composed by P.I. Tchaikowsky in 1891 [55] , 

The results obtained in the present paper show that the "keyboard" of rotating elastic bodies of revolution, 
among which are the glass harmonica and the disc of a brake, is formed by the nodes of the spectral mesh, 
corresponding to angular velocities in the subcritical range. The frictional contact is the source of dissipative 
and non-conservative forces, which make the system unstable in the vicinity of the nodes and force a rotating 
structure to vibrate at a frequency close to the double frequency of the node and at the angular velocity close 
to that of the node. These conclusions agree with the results of recent experiments with the laboratory brake 
[441 1451 l4l)] . The higher modes are more efficiently damped than the lower ones, and a particular frequency is 
selected by the speed of rotation and the loading conditions, including such parameters as the size of the friction 
pads and their placement with respect to the disc. 

It is known that dissipative and non-conservative forces may influence the stability in a non-intuitive manner 
[H [71 [T31 HU [221 EU HH HH EOl E2 ■ We have shown that the former create the latent local sources of instability 
around the nodes of the spectral mesh (bubbles of instability) , while the latter activate these sources by inflating 
and destructing the bubbles. It turns out that the eigenvalues of the damping matrix control the development of 
instability. For better stability both of them should be positive and stay far from zero. If one of the eigenvalues 
of the damping matrix is close to zero or becomes negative, the instability can occur with the weaker non- 
conservative positional forces or even without them. 

With the use of the perturbation theory of multiple eigenvalues we have obtained explicit formulas describing 
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the deformation of the spectral mesh by dissipative and non-conservative perturbations. The trajectories of 
eigenvalues are analytically described and classified. The approximations of the domain of asymptotic stability 
are obtained with the use of the derivatives of the operator and the eigenvectors of the double eigenvalues 
calculated at the nodes of the spectral mesh. Singularities of the stability boundary of a new type were found 
and their role in the development of instability was clarified. 

The theory developed seems to be the first analytical explanation of the basic mechanism of friction-induced 
instabilities in rotating elastic bodies of revolution. The application of the theory to more complicated discrete 
and continuous models of squealing disc brakes, singing wine glasses, and calender rolls will be published 
elsewhere. Comparing the generic results obtained in the two dimensional case with that of the study of the 
rotating string with the singular source of momentum and energy we have shown that the traditionally employed 
concept of a point-wise contact leads to the semi-definite non-generic perturbation operators which suppress 
generic instability mechanism causing squeal. The further progress in the squeal simulation seems to strongly 
depend on creation of the models of distributed contact that do not allow for the semi-definite degeneration. 
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